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TECHNICAL NOTE 3286 


GENERALIZED INDICIA! FORCES ON DEFORMING 

RECTANGULAR WINGS IN SUPERSONIC FLIGHT 

By Harvard Lomax, Frank lyn B. Fuller, 
and Loma Sluder 


SUMMARY 


X 

A method is presented for determining the time-dependent flow over 
a rectangular wing moving with a supersonic forward speed and undergoing 
small vertical distortions expressible as polynomials involving spamd.se 
and chordwis^ distances. The solution for the velocity potential is 
presented in a form analogous to that for steady supersonic flow having 
the familiar "reflected area" concept discovered by Eward. Particular 
attention is paid to indicial-type motions and results are expressed in 
terms of generalized indicial forces. Numerical results for Mach numbers 
equal to 1.1 and 1.2 are given for polynomials of the first and fifth 
degree in the chordwise and spanwise directions, respectively, on a wing 
having an aspect ratio of 4. 


INTRODUCTION 


One of the basic problems arising in the analysis of wing flutter 
boundaries is the calculation of the aerodynamic forces on wings under- 
going small but arbitrary spanwise and chordwise distortions. When the 
wing aspect ratio is large (actually, when the distance between spanwise 
nodal lines is large), these forces are usually estimated by some strip 
theory in which the loading on each spanwise section is approximated 
from that on a two-dimensional wing having the same chordwise distortion. 
This report is concerned with low-aspect -ratio rectangular wings for 
which tip effects are important and the full three-dimensional theory 
must be used. 

The exact linearized solution for the forces on thin rectangular 
wi ngs ( l imi ted, however, to the range where effective aspect ratio 
(a/ M 2 -l A) is > l) travel ing at supersonic speeds has been presented by 
both Gardner (ref. l) and Miles (refs. 2 and 3) in terms of multiple 
integrals involving arbitrary surface undulations . However, the use of 
such solutions in evaluating, numerically say, the forces induced by 
specific wing distortions still presents some difficulties. It is the 
purpose of this report to discuss certain techniques that can simplify 
the labor involved in these calculations and to present numerical tables 
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for the forces induced "by a class of surface deformations, a class gen- 
eral enough to represent the first few mode shapes of rectangular 
plates . 

Mathematic al l y the problem is to find and analyze a solution to the 
four-dimensional wave equation 

toe + V + O’zz - TJ *t't' ‘ 0 (la) 

(where clq is the speed of sound, t* is the time, and x,y,z are space 
coordinates) that satisfies the appropriate boundary conditions. The 
particular form of the solution to he analyzed differs from those pre- 
sented hy Gardner and Miles hut its development is based on the method 
due to Gardner. 

Hadamard (ref. h) studied a generalized form of equation (la) in 
■which the number of dimensions -was arbitrary. His solutions to these 
generalized equations are fundamentally different, depending on whether 
the total number of dimensions is odd or even. In fact, the methods 
Hadamard developed apply directly only to equations for which the total 
number of dimensions is odd. Solutions for the even cases (such as 
eq. (la)) are determined by a ,l method of descent 11 ; that is, the solution 
for the next higher odd-dimensioned equation is found and then reduced 
by (made independent of) one dimension. It is apparent, however, that 
such a technique is in itself by no means unique. Thus, Hadamard found 
the solution to equation (la) by descending from a solution to the 
equation 


+ ^yy + ^zz + ^t't' = 0 

a Q 

but there are many other partial differential equations and groups of 
partial differential equations governing a five -dimensional (x,y,z,£,t) 
space all of which satisfy equation (la) in a plane £ = constant. 
Gardner discovered a set of equations containing equation (la) in a 
f- = constant plane which are simpler than equation (la) in that solu- 
tions could be found and adapted to the boundary conditions for time- 
dependent motion by methods well known to aerodynamicists who have 
studied the flow about wings in steady supersonic flight. This is the 
essential part of Gardner’ s contribution and it represents the tech- 
nique upon which the development of the solution presented in this 
report is based. Actually, Gardner first applied a Lorentz transforma- 
tion to equation (la) and then used his method outlined above. The 
application of such a transformation is unn ecessary and has the dis- 
advantage that the resulting coordinates have lost their direct physical 
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significance. We -will apply Gardner’s method of descent directly to 
equation (la) and then proceed to analyze the solutions so obtained. 

In order to simplify the analysis as much as possible,, we will 
limit solutions to the plane of the wing, and, further, consider only 
indicial-type boundary conditions; in other words , unsteady motions in 
which the wing attains instantaneously, at time zero, a certain span- 
wise and chordwise distortion which is thereafter fixed. It is well 
known that the transient responses to these indicial motions can be used, 
in a superposition integral, to obtain responses to many other types of 
unsteady motion; in particular, responses to the harmonic oscillations 
of nonrigid wings. 

Finally, the principal interpretation of the results will be made 
in terms of generalized forces, since these can be used directly in 
either flutter or gust studies, and it will be shown that the amount of 
labor required to calculate such forces is reduced by using reciprocity 
relations derived from the general theorems presented in reference 5. 


LIST OF IMPORTANT SYMBOLS 


A aspect ratio 

ao speed of sound 

&i n amplitude of indicial-downwash distribution (See eq. (2a).) 

B(p,q) beta function (See eq. (B15a).) 

B x _ x 2(p,q) incomplete beta function (See eq. (B15b).) 

C(x!,y!) influence function for effect of side edge (See eq. (A10).) 

lift 


C L 


lift coefficient. 


<lo s 


'La 


indicial lift coefficient due to angle -of -attack change, 

SCt 

without pitching, Ct = — — 

oa 1 ot=0 

indicial lift coefficient due to pitching for a wing rotating 

Sc T 

about its leading edge, Cr. = -r-= 

q oq 


q=0 


'm 


pitching -moment coefficient, positive when trailing edge 
tends to sink relative to leading edge, - Qgie - ^ 
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indicial pitching-moment coefficient due to angle -of -attack 
change (without pitching) measured about the leading edge, 

n _ ^Cm 

a “ -s 

da a =o 

S’ 

indicial pitching-moment coefficient due to pitching measured 
about the leading edge for a wing rotating about its lead- 


ing edge, Cjn = 

^ oq 

q=0 

c 

wing chord 



generalized indicial force coefficient (See eq. (36).) 


generalized indicial force coefficient (See eq. (37)0 

h(x,y,t) 

distance of wing camber line from z = 0 plane 

M 

Mach number 


Ap 

loading coefficient (pressure on the lower surface minus 
pressure on the upper surface divided by free-stream 
dynamic pressure) 

CD 

binominal coefficient. 

(A = 

\ m / ml (n-m) 1 

q 

• 

c Q 

dimensionless rate of pitching, — 

Uo 

*o 

1 2 

free-stream dynamic pressure, — P 0 U 0 


generalized coordinate 


Q r 

generalized force corresponding to the generalized coor- 
dinate q r 

R.P. 

real part of 


*o 

./(x-x-l) 2 + (y-y x ) 2 


r i 

J (x-x 1 ) 2 + (y+y x ) 2 



*c 


J (x-x 1 ) 2 -p 2 (y-y ;L ) 2 



MCA TN 3286 


5 


S 

S a 

•S c 

t 

t' 


"tm 

T 

U 

Uo 

W 


w 

x,7,z 


X 3>y 3 } ^3 
x 4 ^y 4 ,t 4 


wing semi span 
wing area 

area of acoustic plan form 

area of reflected acoustic plan form 

a 0 t' 

time 

t 

c 

x + Mt 

P 

wing kinetic energy • 
wing potential energy 

forward speed of wing 

(St) 

\ 8z / 2=0 

vertical velocity 

Cartesian coordinates, fixed relative to the fluid at 
infinity 

coordinates with origin on center of wing leading edge 
(See sketch (2).) 

coordinates with origin on center of wing leading edge at 
time zero (See sketch (n).) 




K-m 


X x (t|) 


X 

c 

Mx + t 

P 

- ytm 2 - t) 2 ^ 
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a angle of attack (angle "between flight path and plane of 

wing), radians 

3 a/m 2 - 1 

9 wing angle of pitch relative to horizontal, positive when 

trailing edge lies "below leading edge, radians 

| coordinate measuring fifth dimension 

P Q free-stream density 

cp velocity potential 


9 


(1) 


portion of velocity potential induced *by sources in acoustic 
plan form 


<P 


(2) 


portion of velocity potential induced "by presence of side 
edge 


t 


potential function in five -dimensional space- 

to 


Subscripts 


A,B,C regions in an x,£ plane (See sketch (d).) 

u upper side of wing, z = 0+ 

1 singularity (e.g., source) position 

I,II,...VIII regions on wing shown in figure 1 


STATEMENT OF THE PROBLEM 
The Governing Equation 


Assuming a wing’s vertical motion is of such a nature that the 
velocities induced in the fluid are small relative to the magnitude of 
the wing’s steady forward motion, the normalized form of equation (la) 


cp + 
xx 


cp + cp 


zz 


- ftt - 0 


(lc) 
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where t = a^t', can be used as the governing partial differential equa- 
tion of the flow field. This equation applies to the determination of 
the velocity potential when the body or wing in question moves through 
the fluid, the axes remaining fixed with respect to the still fluid 
infinitely distant from the origin. For convenience we place the wing 
leading edge on the y axis at t = 0 and the side edge on the x 

axis. The wing flies at a constant f orward (in the negative x direc- 
tion) speed so at subsequent times the leading edge lies along the line 

x = -Mt, where M is the Mach number, and the side edge moves along 

the x axis as shown in sketch (a) . 

The Boundary Conditions 


The fluid velocity normal to the 
surface of a solid moving in a friction 
less fluid must be zero. If the equa- 
tion of the solid's surface is repre- 
sented by 

G(x,y,z,t* ) = 0 

this boundary condition can he expressed Sketch (a) 

mathematically, in terms of the coordinate system used in equation (lc)> 
as 

bt 1 bx dx by by bz bz 

Consider a thin surface near the z = 0 plane- The equation of the 
camber line of this surface can then be expressed in the form 

G(x,y,z,t’) = z - h(x,y,t' ) = 0 

and, assuming that thickness and lifting effects can be separated lin- 
early, the boundary condition for the camber line becomes 

8h dq> 8h . cicp dh _ cKP _ n 

St 1 5x c)x by by bz 

If the derivatives of h -with respect to each of the coordinates are 
small, the two middle terms can be neglected and the expression for the 
boundary condition reduces to 

— =w u (x,y,t») 

St ' Sz z =0 
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We wish to simulate a rectangular wing deformed indicially hy tend- 
ing in the spanwise and chordwise directions. For this purpose, on the 
portion of the z = 0 plane occupied hy the wing plan form, the vertical 
velocity, which determines the wing shape according to the previous 


t < o 


t > o 


where c is chord length, a^ n is a constant and l and n are 
integers >0. 

The expression (x + Mt) ^ is used so that for l > 0 the tangent 
to the wing camber line at the leading edge is tangent to the flight-path 
angle of the leading edge. Consider, for example, the case l = 1, n = 0. 
The downwash 


equation, is assumed to have the form 


w u = < 


XX-M^T 


w u = (x + Mt) 


represents an infinite class of surface shapes having the form 


h(x,y,t) = —12- [(x + Mt) 2 + f (x,y) ] ( 2 ) 

2cU 0 v 

where f (x,y) is an arbitrary function and h is, hy definition, the 
distance of the wing’s camber line from the z = 0 plane. Since, within 
the accuracy of linearized theory, the solution for the flow about the 
wing depends only upon the value of v u (x,y,t) , the loading on all the 
wings represented hy the above equation is the same. 

Let us inspect the tiro special cases 

(i) f(x,y) = -x 2 
(ii) f (x,y) = 0 


For case (i) 


h(x,y,t) 


a 10 M 

2cU 0 


(2xt + Mt 2 ) 


u 
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and the wing is a flat plate pitching at a uniform rate about its lead- 
ing edge -which is following the flight path 



a 10 M 2 t 2 

2cU 0 


as shown 1 in sketch (h) . Hence, at time 
path of the leading edge is 


t the tangent to the flight 



d ( h W dt * a io t 


-U r 


The slope of the leading edge of the plate at the same time is 


Was 

and the two slopes are seen to he equivalent. 
For case (ii) 


h(x,y,t) = — — ( x + Mt) 2 
2cU 0 


1 The z scale in both sketches (b) and (c) is purposely distorted 

in order to make the drawings clear. A basic assumption used in setting 
up the boundary-value problem, by means of which the loading was deter- 
mined, was that the surface of the wing must remain near the z = 0 plane . 
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and the wing is a plate which obtained a sudden parabolic camber at 
t = 0, a shape it maintained thereafter as shown 2 in sketch (c). 


t - 

1 

l 

0 

TTtr^ 

» z 

f, 

rtrw 

- J 

rrw 1 

iTrrr^^ distribution 

X 

^ 1 





flpUL -q a f leading edge 

dt 

Sketch (c) 


The problem is linear, so it vill he sufficient to determine a 
solution for arbitrary l and n, and then add results for any combina- 
tion of terms as desired. Thus, the complete boundary conditions to 
be studied are 


w u (x,y,t) 



(2a) 


over the wing plan form, and, since the loading is zero over the remain- 
ing portion of the plane 


dcp 

St 


0 off the wing 


z=0 


since the loading is given by 


Ap = _ k _ z'&p'n 
% ~ UqM V3t J z =o+ 


S0LUTI0N FOR THE POTENTIAL 


(2b) 


Figure 1 shows the wing plan form on the surface of which the 
potential is required, together with the system of axes; also, traces 
in the z = 0 plane of the wave system set up by the indicial motion of 
the wing are indicated. The wave pattern for only tiro edges is shown; 


2 See footnote 1 on p. 9* 
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the flight speed is supersonic so the trailing edge has no effect on 
the velocities induced over the wing surface, and the results are valid 
(in their entirety) only for {3A > 1, so the opposite edge either has 
no effect or one that can he incorporated hy simple superposition. 

The wave traces divide the wing area into several regions, indicated 
hy the Roman numerals, in each of which the analytical formulation for 
the potential is different. Region I consists of that part of the wing 
where the effect of neither the side edge nor leading edge has yet heen 
felt. In region II, the side-edge influence is acting (the line y = t 
is the trace of the starting cylindrical wave from the side edge y = 0) 
hut not the leading edge. Region III is the part within the starting 
cylindrical wave from the leading edge, hut outside the influence of the 
side edge. This region, and region V, are further subdivided for reasons 
that will appear later. Region IV is a compound region; potential there 
can he found hy adding the potentials for regions II and III and sub- 
tracting the potential for region I. Region V consists of the portion 
of the wing within the spherical wave originating at the wing corner. 

The flow over the part of the wing comprising regions VI and VII has 
reached a steady state relative to a point on the wing, and the poten- 
tial there is just that for the corresponding parts of a rectangular 
ling with the proper downwash distribution in steady motion. Finally, 
region VIII is again a composite region, its potential being the sum of 
potentials for regions III and VII less the potential for region VI. 

All the regions just listed, with the exception of region V, are 
actually governed hy the three- (total) dimensional wave equation and 
the potential therein could he obtained hy methods applicable to this 
simpler equation. However, in this report we shall present a -unified 
approach and the problem will he solved hy the same method in all 
regions . 


Review of Kirchhof f 1 s Formula 


The solutions developed in the subsequent sections are more clearly 
interpretable if they are compared with certain known results that have 
already heen determined for the indieial motion of nonlifting wings with 
symmetrical thickness distributions or lifting surfaces with all super- 
sonic edges. The purpose of this section is simply to review briefly 
some of these latter results. 

As in steady-state wing theory, there is a formula for time- ' 
dependent flows that relates the velocity potential to a distribution 
of time -dependent sources and doublets over a certain region in the 
wing plane. This formula is due to Kirchhof f, and some of its aero- 
dynamic uses are discussed in reference 6. Kirchhoff's result is 
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immediately applicable in the study of unsteady lifting-surface problems 
vben the potential can be represented by sources alone, that is, when 
the upper and lower surfaces of the -wing do not interact, as is the case 
in regions I, HI, and VI of figure 1. 

Kirchhoff’s formula for source distributions can be written 

■p(x,y,0,t) = - i dy x (3) 

Sa 

where 

r Q 2 = (x - xj 2 + (y - yj 2 


The brackets on indicate that the retarded value is to be taken 


KJ = VuU^y^t-rc,) 


and S a indicates that the region of integration is the acoustic plan 
form corresponding to the event (x,y,0,t) . These concepts are discussed 
at length in reference 6. 

As has been pointed out, equation ( 3 ) holds for each of the regions 
I, III, and VI, but the area of integration S a differs considerably 
from one of these regions to another. Consider, for example, the deter- 
mination of cp for region III, denoted cp^ ^ . Part of the boundary of 
the acoustic plan form Sa, is found by eliminating T between the 
equation of the leading edge, x x = -MT, and the expression 

(x - xj 2 + (y - yj 2 = (t - T) 2 


which gives the outer boundary, at "time" t, of all the disturbances 
that, operating at "time” T, can produce an effect at the point (x,y) . 
This boundary is the ellipse 


where 




(y - yj 2 = V 2 


„ _ Mx+t x+Mt 

% “ o ; % - 


(4a) 


0 


0 
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If tlae point (x,y) lies within the cylindrical wave from the leading 
edge, that is, -t < x < t, the ellipse of equation (4a) comprises only 
part of the acoustic plan form, the remainder being hounded by so much 
of the circle 


(x - xj 2 + (y - yj 2 = t 2 


(4b) 


as lies on the wing at time zero. Sketch (d) shows the three possible 
acoustic plan forms for points in region III. The limits for the three 
types are 


(i) t > x > 0 

(ii) 0 > x>-t/M 

(iii) -t/M > x > -t 



and these correspond to the subregions III a , IHb, and IHc identified 
in figure 1. Using equation (3)> we can write the potential in, say. 



l4 


KACA TN 3286 


region IH a as 

y+t x+V t 2 -(y-y 1 ) 2 

**■-= / M ^ 

y-t x-Vt 2 -(y - 7l ) 2 


9- 


dxi + 


_1_ 
2 jt 


y+*/ t 2 -x 2 X 1 (y-y 1 ) 

/ / 

y-*/ t 2 -x 2 x— a/ t 2 -(y-y x ) ; 


KiJ 

=• r O 


dx. 


(5) 


where 


L (y - Xi) = ^ --/tm 2 - (y-y x ) 2 J 


Gardner's Method of Descent 


Equation (l) governs a four-dimensional x,y, z,t space. Our object, 
of course, is to find for this equation a solution that satisfies the 
boundary conditions in the z = 0 plane as specified in equations (2a) 
and (2b) . Obviously, we can always construct a space of more dimensions 
governed in an arbitrary way except that it must satisfy equation (l) in 
an x,y,z,t hyperplane. Then, if a solution in this higher dimensional 
space which satisfies equations (2a) and (2b) in the x,y,z,t plane can 
be found, it represents for £ (the additional dimension) equal to some 
constant the solution to our problem. This characterizes the method of 
descent. It is not obvious, of course, that such a method leads to any 
simplification; but, with a proper choice of the governing equation for 
the new space, such a possibility always exists. 

There are examples where various applications of this method have 
proved to be useful. Hadamard* s use of the method, mentioned in the 
introduction, is classical. A simple application of his method is the 
derivation of the velocity potential for a source in a two-dimensional 
supersonic flow field. This potential field (which amounts to a step 
function, the step occurring at the Mach wave) is easy to derive if one 
considers a three-dimensional field with a line of sources normal to the 
free stream and uniform in strength. The two-dimensional field mentioned 
above follows Immediately by descent. 

In other examples the additional dimension is measured with imagi- 
nary numbers and the additional law for the extended space is the 
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requirement that the functional dependence on the resulting complex 
variable shall he analytic. The method of descending in the latter 
case is associated with the study of analytic continuation. In partic- 
ular, Riesz's method (discussed in ref. 7) for solving equation (l) 
illustrates these concepts. 

Gardner's method for solving equation (l) is to define a five- 
dimensional space in which a potential function i|r is governed by the 
. equations 


+tt " *xx - = 0 

♦|| - V “ 0 (6b) 


and show that solutions to equations (6) in this space are general 
enough to contain general solutions to equation (l) in a plane 
£ = constant. We shall, therefore, proceed by analyzing these equa- 
tions and eventually let £ approach a plane in which the boundary 
conditions of equations (2a) and (2b) are satisfied. For convenience, 
the latter plane is taken to be the £ = 0 plane. 

Since equations (6a) and (6b) are linear, a number of possibili- 
ties exist for the choice of the dependent variable 4 r (x,y,z,0,t) . 

Aside from the more obvious choice ijr(x,y,z,0,t) = cp(x,y,z,t), where 
is the velocity potential of equation (l); for example, one could let 
ijf(x,y,z,0,t) = cp x (x,y,z,t) or again, t^(x,y,z,0,t) = cp(x,y,z,t) . These 
various choices amount only to relatively minor differences in the 
detailed technique of the subsequent analysis. If, in imposing the 
boundary conditions of equations (2), one is to use only source-type 
solutions for both equations (6a) and (6b), the last choice is suffi- 
cient. Therefore, set 


~_ 8 _ 

81 


i(x,y,z,|,t) 


J £=0 


<p(x,y,z,t) 


(7) 


Wow differentiate equation (6a) with respect to z and set 3 z = 0. 

3 

It can be shown that the solution satisfies the equation 


lim f lim 
z— >0 \_£— >0 



(x,y, z,|,t) 


lim 
z — >0 


<p(x,y,z,t) 



tg(x,y,z,£,t) 


} 
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Defining 


W(S,x,y,t) 

oz z =0 


(8) 


equation (6a) can be expressed in the form 


w tt " w xx " % - 0 


(9) 


and the boundary conditions for equation ( 9 ) are given directly by equa- 
tions (2) . Thus on the -wing 


6w 

_ 89 


51 

1=0 

z=0 

and off the wing 




cw 



St 


= v u (x,y,t) = a 


In 


^x+MtV 


)(* 


(10a) 


1=0 


= cp t (x,y,0,t) = 0 


(10b) 


Assuming equation ( 9 ) to have been solved for the boundary condi- 
tions given, by equations (10), we return to the second of the set of 
partial differential equations (6) , specifically. 


♦st - V - - 0 


Prom equation ( 8 ), it is seen that the solution to equation ( 9 ) 
yields the result 


5pjr 

5z 


z=0 


= known function of 


7>i 


on the wing 


Further, the boun dar y conditions for the original 'problem in (x,y,z,£,t) 
space require that 9 be an odd function with respect to z, and con- 
tinuous across the z = 0 plane except over the wing plan form. Thus 9 
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must tie zero for z = 0 except over the wing plan form. The continuation 
of this condition into (x,y,z, 5,t) space then implies, according to equa- 
tion (7)> that off the wing 


at 

85 


l z=0 



Hence, hoth the second partial differential equation and its houndary 
conditions are identical in form to the first set given by equations ( 9 ) 
and (10), respectively. Applying equation ( 7 ) to their dual solution, 
we obtain the desired result 


_85 


t(x,y,0,5,t) 


J5=0 


<P(x,y,0,t) 


for the potential on a rectangular wing (with PA > l) in supersonic 
unsteady motion. 


The General Expression for the Potential 


The method outlined in the preceding section will now be applied 
to obtain integral expressions for the potential, in any region of the 
rectangular wing shown in figure 1. Consider first equation ( 9 ) for 
W(£,x,t). This equation is the same partial differential equation as 
that which governs supersonic steady flow. Further, the boundary values 


in the £,x,t space are identical to 
those representing a thin planar wing 
in a steady supersonic flow. Since the 
Mach number in the steady-flow analog 
is */~ 2 , the equivalent plan form of 
this wing (shown in sketch (e)) is a 
sweptf orward wing tip having all super- 
sonic edges (i.e., the component of the 
free-stream velocity norma], to all edges 
is supersonic) . 

Since all edges of the equivalent 
wing plan form are supersonic, the 
solution for W can be written imme- 
diately in terms of "sources" only, 
their strength being given by equa- 
tion (10a). Thus, by analogy with 
the well-known results of supersonic 



Sketch (e) 
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wing theory, -we have 


W(i,x,t) 



w^Xj+Mt^y) dxidti 

'/(t-t 1 ) 2 -g 2 -(x-x 1 ) 2 


( 11 ) 


where r is the area on the wing cut out hy the forecone from the 
point (i,x,t), see sketch (e). The analytic form of W will differ 
considerably in each of the three regions above the equivalent wing 
shown in sketch (f ) . 


The value of _ W given by 
equation (ll) now becomes a bound- 
ary condition for the solution of 
equation (6b). Thus, over the 
portion of the z = 0 plane for which 
y > 0, | > 0, the variation of 

8*1 



8z 


is now known and for 

'z=0 

y < 0, | > 0 
8 * 


the condition 


81 


= 0 applies. (These condi - 


z=0 


stream. 


Sketch (f) 

As shown in sketch (f ) , 


tions are still not sufficient to 
determine a -unique solution unless 
the further restriction is imposed 
that the loading falls to zero as 
the edge y = 0 is approached, 
i.e., as y— >0 +.) Again we observe 
that these boundary conditions and 
the partial differential equa- 
tion (6b) are identical to those 
studied in connection with a sta- 
tionary planar wing in a supersonic 
solutions from the t,x,| space above 
as W A , W B , and W c , depending on the 
t = constant plane. Sketch (g) shows 


the |=0 plane are referred to 
relation between x and g in a 

the five different boundary-value problems formed by the various com- 
binations of W A , ¥3, and W q occurring along constant x lines in the 
x,S plane; and the corresponding regions in figure 1 for which each 
applies. Each of these five problems is directly a n alogous to the 
boundary-value problem encountered in steady-state lifting-surface 
theory, of a planar, rectangular lifting surface in a steady supersonic 
stream. The "leading edges" of the se anal ogous rectang ular plan forms 
lie along the lines = t, = J t 2 -x 2 or depending on the 

value of x, and the "side edge" lies along the line y = 0. Hence, by 
means of this steady-flow analog, we can immediately write the solution 
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\ 

(Hi) -t/M*x*Oi M h JT b 


to equation (6b) in the form 

^,7,0, |^t) = 


19 


W(x,y 1 ,i 1 ,t)dl 1 dy 1 

- (y ^) 2 


•where only the area of integration 0 
must he discussed. 

Two possibilities exist for the 
shape of 0 . First, if the point 
I, y lies to the right of the dashed 
lines in sketch (g), which in the 
analogous steady-flow problem repre- 
sent the traces of the Mach cones 
from the leading-edge tips, 0 is the 
triangular area shown (for region 
Ilia) in sketch (h) part (i). If 
however, £, J lies between this line 
and the side edge, y = 0, 0 is the 
trapezoidal area shown (for region V a ) 
in sketch (h) part (ii). The latter 
is a well-known result used in steady 
supersonic lifting-surface theory and 



(ii) 0<x<t , y< /t e -x 

Sketch (h) 
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first developed "by Evrard (ref. 8). The division of the five kinds of 
problems illustrated in sketch (g) into the final twelve, represented by 
the regions in figure 1, is brought about by the various combinations of 
W^, Wjg, and ¥q that can occur in the area a as the point £, y assumes 
all necessary values on the wing. 

When ijr has been determined, the potential in the physical plane is 
found by equation (7)> or, combining equations (ll) and (12), 


cp(x,y,0,t) = 


1_ ^lim pp pp ^u ( X i + Mt 1 , y j_) dXj.lt x 

t 2 85 JJ J (|_| i )2_( y _y i )2 JJ J (t-t 1 ) 2 -i 1 2 -(x-X 1 ) 2 


(13) 


A detailed analysis of equation (13) for a point x,y,t in region 
V a of figure 1 is given in Appendix A, and a study of this analysis 
enables one to write the results for all regions without difficulty. > 


Interpretation of the Results 


The results of the rather involved analysis given in Appendix A 
can he interpreted in terms of the known solutions for simpler boundary 
conditions. These latter solutions have already been reviewed in a 
previous section in which it was shown that the potential on a lifting 
surface with all supersonic edges can be written in the form 

■ - * ss*** 

Sa 

From Appendix A it is found that the potential at a point on a rectangu- 
lar lifting surface can always be expressed as the sum of two parts 


cp(x,y,0,t) = cp^ 1 ^ (x,y,0,t) - (x,y,0,t) (l4) 

where 

Sa 

and 

q/ 2 \x,y,0,t) = - i ff 

s c 


C(x 1 ,y 1 )dx 1 dy 1 


(15b) 
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The value of Cfe^yj.) is given by equation (A10) in Appendix A and the 
areas of integration, S a and S c , are illustrated for the various regions 
I through VIII in figure 2. 

Let us first inspect equations (15) in light of their possible 
analogy with the familiar solution for the steady-state, rectangular 
lifting surface. If a rectangular wing having arbitrary twist and camber 
is placed in a steady supersonic flow, the solution for the potential on 
its surface can also be expressed as the sum of two parts 

<p(x,y,0) = q)(^(x,y,0) - <P^ (x,y,0) ( 16 ) 


where , if 

r c 2 = (x - x^ 2 - p 2 (y - yj 2 


and 


cp^(x,y,0) 


1 rr 

* JJ r c 

Sr 


(17a) 


<p^(x,y,o) 


1 rr w-qdXj.dy 

ft <J J Tq 

S 2 


(17b) 


These equations can be construed in 
the following simple way: Equation 

(17a) represents the potential induced 
at x,y,0 by a distribution of 
sources over the wing plan form, each 
source having a strength proportional 
to the local streamwise slope of the 
upper surface. The area Sj_, as 
shown in sketch (i), is the portion 
of the wing within the Mach forecone 
from x,y,0. Equation (17b) has a 
similar interpretation; it also rep- 
resents a distribution of sources 
over the wing, each having a strength 
proportional to the local slope of 
the upper surface. But the area of 
integration S 2 is now that portion 
of the wing within the Mach forecone 
from the point x, -y,0; that is, 
within the cone which forms a mirror 
image of the physical Mach forecone 
in the vertical plane containing th.e 





u 



Sketch (i) 



22 


HACA TN 3286 


"wing’s side edge. The potential 9 2 (x,y,0) represents the difference 
"between the potentials for a wing with a vertically symmetrical thick- 
ness distribution and a surface with no thickness having the same shape 
as the upper surface of the nonlifting wing. 

Let us return now to equations (15) . Just as in the steady-state 
case, cp^ 1 ) (x,y,0,t) represents the potential induced at x,y,0 "by a 
distribution of sources (see eq. (3)) over the wing plan form, each 
proportional to the local slope of the wing, but now, since the wing is 
in motion, with the added condition that they be local slopes at the 



Sketch (j) Sketch (k) 


appropriate time. The area S a , shown in sketch (j), is just the 
acoustic plan form defined earlier in the discussion of equations (3) 
and (^) . Physically, Sa represents those points on the wing from 
■which disturbances can, at the time t, influence the flow at x,y,0. 

It is the generalization, in the stationary coordinate system, of the 
wing area bounded by the Mach forecone. 

The relation between cp^ 1 ^ (x,y,0,t) and cp^ 2 ^ (x,y,0,t) is similar 

to that between their steady-state analogs. Thus, again, cp^ 2 ) (x,y,0,t) 
represents the difference between the potentials for an uncambered non- 
lifting wing and a lifting surface having the same shape as the top of 
the nonlifting wing. A more striking similarity lies in the relation 
between S a and S c . 

We have already seen that Sa, is the acoustic plan form, and, as 
it turns out, S c jLs the reflection of the acoustic plan form (see 
sketch (k) ) in~he vertical plane containing the side edge - a situation 
identical to that existing between S x and S 2 in the steady-state case. 
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(In other words, Sa is the acoustic plan form for the event x,y,0,t, 
and Sc is the acoustic plan form for the event x, -y,0,t.) Physi- 
cally, S c represents the portion of the wing’s lower surface contain- 
ing disturbances which can, at the time t, influence the flow at 
x, y,0 on the wing’s upper surface. At this point the similarity 
"between the steady and unsteady solutions ends since the influence of 
the slopes in the reflected plan form is not the same as it is for the 
slopes in the basic acoustic plan form; the influence in the former 
case now being given by the integral C(x 1 ,y 1 ) defined in equation (A10) 

One can show, by simply referring the results given in equations 
(15) to a coordinate system fixed on the wing, that equations (15a) and 
(15b) are identical, respectively, to equations (17&) and (17b) -when 
they apply to regions VII -and VT in figure 1; regions in which, for 
indicial-type motions, the flow is steady relative to the wing. Hence, 
equations (15a.) and (15b) extend Eward’s n reflected area" concept to 
all parts of a rectangular wing in supersonic unsteady motion. 4 


TTTP GENERALIZED FORCES 
Review of Lagrange’s Equations of Motion 


In order to define more clearly the subsequent concepts and nota- 
tion, we will briefly review Lagrange’s equations of motion as applied 
to distorting wings and will examine a simple application to a rectan- 
gular wing. 

Lagrange's equations are usually written 


d dT _ _8T_ _dU 
dt’ dq r dq r Sq r 


Q r ; r = 1,2, 


(18) 


where 

T kinetic energy of the wing 
U potential energy of wing 
Q r a generalized (external) force 
q r a generalized coordinate 


4 It is of further interest to notice that equation (15b) can be 
reduced to a double integral involving w u (^,y x ) by using, for example, 

the transformations £ = x x + Mt x and t = t - 1, and integrating in the t 
plane . 
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In the present application q r is the amplitude at a given time of a 
polynomial measuring h, the. vertical displacement of the •wing' s camber 
line from the z = 0 plane. Thus, relative to an x 3 , y 3 coordinate 
system that is fixed on the wing, see sketch (Z) 


Mx 3 ,y 3 ,t*) 


^ q s (t«)P s (x 3 ,y 3 ) 
s 


The wing 1 s kinetic energy can "be "written 


(19) 


T = JT i h 2 m(x 3 ,y 3 )dx 3 dy 3 (20) 

S 

where m is the wing mass per unit plan-form area. Using equation (19)* 
we find 


d 8 t 

dt T 8q r 


l 



Pr( x 3^y 3 )Ps( x 3^y3)m(x 3 ,y 3 )dx 3 dy3 



= 0 


} ( 21 ) 


The potential energy is usually 



difficult to evaluate analytically. 
However, it can often he determined 
experimentally (as will he seen) hy 
measuring the frequencies of the free 
vibration modes.. For the present 
assume that the wing is a homogene- 
ous plate of constant thickness. 

The potential energy for such a 
wing can he expressed as (ref. 9) 



Sketch (Z) 
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which leads to the equation 


= D Z % f f V2PlV2Ps " 2(1 " ( 


w 1 S 2 P r S 2 P s 

2(1 - d) [ ± r 8 


2 ^ 


l 3 2 p r 8 2 p s 8 2 P r £ 2 P e 


^3^3 


2 dx 3 2 by 3 2 8x3^7 3 Sx 3 £y 3 ^ _j 
where d is Poisson’s ratio, V 2 = 8 2 /8x 3 2 + 5 2 /Sy 3 2 , and 


(23) 


D = 


2 (Young’s modulus) (plate thickness) 
3(1 - d 2 ) 


Now, if the generalized coordinates have "been normalized so that 
each measures the amplitude of a free vibration mode, all terms in equa- 
tions (2l) and (23) involving the integral of the product of P r and P s 
are zero. Assuming, henceforth, such normalization , we can write 


<3 r 



P r 2 ( x 3^Y3) m ( x 3^y3) dx 3 d y3 + D( lr 



2(1 -d) 


" 8 2 P r 8 2 P r 
-dx 3 2 8y 3 2 


( 


£ 2 P r 

dx Q dy 3 



dx 3 dy 3 = 


Qr* > 


r = 1, 2, 


(2*0 


Finally, dividing through by the coefficient of q r and expressing a 
generalized force as the integral, over the wing plan form of the product 
of the rth mode shape and the loadings 5 s(Ap) s induced on the wing by 
each of the mode shapes considered, we find 


q r + q r ^ 2 



// IV 2 (x 3 ,y 3 )m(x3,y3)dx3dy 3 
S 


where io r is the frequency of the rth' free vibration mode . 


(25) 


will write (Ap)_ = q 0 (Ap/q 0 ) g where q,-, is the free-stream 
dynamic pressure. This Is possible without a confusion of notation since 
the generalized coordinates are expressed as q 1 ,q 2 ,q 3 , • . • and 
exclude the term q^ 
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If the free-mode frequencies are experiment ally determined, equa- 
tions - such as equation (23) - giving the wing* s potential energy, 
never have to he evaluated. Further, in such cases, equation ( 25 ) 
applies to quite general wing structures with varying density. Usually 
in the application of equation (25), one uses the actual frequency w r 
of the free mode hut, in evaluating the aerodynamic forces, uses an 
analytical expression that only approximates the rth mode shape. Let 
us examine the generalized force term in equation (25), taking, for 
simplicity, only one term of the sum; 


Or = 9° J'J M*3>y 3 ) (j^) dx 3 a y3 


( 26 ) 


According to what has gone before, the mode shape polynomial P r (x 3 ,y 3 ) 
has the form 


.(w a )=(^y(^) g 


(27) 


•while (Ap/q 0 ) s is the loading coefficient corresponding to an indicial 
deflection (see previous section on boundary conditions) 


* - £ ,3(1) 


*) B [( 




\l+i 


f 


X, - Mto 


(28) 


•which gives a vertical velocity distribution 


w u = u o ^sd) ( -fj 


(29) 


Now a generalized indicial force coefficient can be defined as follows: 




dXgdyg 


(30) 


> in / 


(The calculation of these quantities f?“(t*) will be elaborated in the 

JS 

next section.) Since the generalized force Q r is intended to apply 
to any motion, not necessarily indicial, it is necessary to apply 

7 n 

Duhamel’s integral to the indicial force coefficient f^ (t’)j thus. 


is 


a n v 


q s (D 


dr* 


( 31 ) 
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As an example , consider nov the simple one degree of freedom 
vibrating plate illustrated in sketch (m) . The plate is fixed to the 



val.1 and restrained along its leading edge. The mode shape is assumed 
to have the form 


/ \ 2 / v 

h = q (t ») ( fa ) f Za ) 
c J \ C J 


(32) 


bo for a plate with uniform density and thickness 


msc / ^ 

25 V c 


m J 6X3 ^ 2 ( x 3^ y s) 

-s o 

Equation (25) now becomes 

«X ♦ <o 1 2 q 1 = -22-(s)\ 

x • msc \c J y 

For this case, we have the generalized indicial force coefficient 
fg^t 1 ), and so 


(33) 


Q i = ^ sc ) f q i (t - T,) 


*£<*’> 


q x (D 


dr* 


(34) 


Therefore, equation (33) can he written 

2 5qn / 


-p 

( 

f 12 (T») 

) ^ y 

O 

22 ' 

q x (l) 


dr« 


(35) 
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The Generalized Indieial Force Coefficient 


It is clear from the previous section that a study of the dynamic 
"behavior of rectangu lar "wings moving at supersonic speeds can "be carried 
out if one can obtain values of the generalized force coefficient, 

fjg(t*), as defined "by equation (30). We will now show how these values 


can "be obtained from the solution to the aerodynamic bo undar y- value 
problem represented by equation (14) . 


tpo 


- x 4 +Mti 

h - U 



Sketch (n) 


It was convenient in developing 
equation (l4) to use a coordinate 
system - x,y, z,t - which was fixed 
in space so that the left edge of the 
wing moved along the x axis as 
shown in sketch (a) . On the other 
hand, in studying the dynamic problem 
it was more convenient to use an 
x 4 ,y 4 ,z 4 ,t 4 system which was fixed 
in space so that the wing 1 s spanwise 
center line moved along the x 4 axis, 
see sketch (n) . Let us first con- 
sider the problem of transferring the 
results in terms of the x,y, z,t 
coordinates to the x 4 ,y 4 ,z 4 ,t 4 
system. 

The indieial force coefficient 

F^ n (t’) is defined as follows: 

J6 


c-Mt 


? 2n 

Jg 


-Mt o 


.g 


«,) (0 (a 


^ln 


(36) 


In order to transfer the axes from the set shown in sketch (a) to the 
more convenient set of sketch (n), so that mode shapes are symmetric or 
asymmetric about the wing's spanwise center line and the force coeffi- 

7 Tl 

cients denoted fj can be determined, we proceed as follows. First, 
the loading coefficient for a wing in the (x,y) system with downwash 




MCA TN 3286 


29 


given by 



H=o 


is obtained. This loading coefficient can be -written as a sum: 



This last integral can be written as 



O |Vj 
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By using equation ( 36 ) 


we find 



X-1 ) eV 

2 


I 


V=o 



n 


p=o 


n-ti 



(37) 


where all forces are responses to a unit indicial disturbance. Note 
that if equation (37) is applied in the case of a wing cantilevered on 
a wall, hoth n and g must "be even in order to satisfy the boundary 
conditions of reflection in the wall. 


By superimposing boundary conditions and their resulting solutions, 

one can further show that the value of f ^ given by equation (37) is 

valid for al 1 reduced aspect ratios 0A greater than 1 in spite qf the 

fact that the value of F^ n given by equation ( 36 ), as it stands, 

J6 

applies only to wings for which PA is greater than 2. 


7 

Given f (t ' ) , one can determine the generalized force associated 

do 

with the generalized coordinate q r by means of the superposition inte- 
gral as illustrated by equation ( 3*0 . 


Details of Calculation 


The details of actually evaluating the indicial force coefficients 
from the solution for the potential presented in the first part of this 
report are discussed in Appendix B. Considerable labor is involved in 
such calculations, and an attempt was made to discover recursion formulas 
by means of which certain derivatives, for the rectangular wing, could be 
expressed as combinations of others. This attempt was successful and 
yielded the following results 

Consider equation ( 36 ) . Integrate the x integral in this equation 
by pails, setting 

u(x) = r B y£ — — dy; dv(x) = (x+Mt)^ dx 

J 0 % 

Then, since by equation (BT) iu Appendix B 

.SL Ap Zn _ l Ap^~ x,n ^ 
cix c qQ 


l > 0 
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one finds 


In _ _Z_ f l-i,n 

is j+il° s j+i^gj 


( 38 a) 


Inspection of equation (37) shows that the same relation holds for the 

2n 

generalized indicial force coefficients f that is, 


Zn Z f_Z-i,n „Z-i,n"1 
is ~ j+1 \ r o g " j+i,g J 


(38b) 


From this relation, it is 

determined hy integration; the 
can he found hy combination of 
shape index i . 


seen that only the forces F 


on 

is 


forces for higher values of the 
results for different values of 


need he 

index Z 
the mode 


As a simple illustration of the results presented so far, we can 
calculate the indicial force derivative for the cases Z = n = g = 0, 
j ss 0, 1. The case j = 0 corresponds to the indicial lift coefficient 
for a flat, sinking, rectangular wing, and the case for j = 1 corre- 
sponds to the indicial pitching -moment coefficient for the same wing. 
Since n = g = 0, equation (37) gives 


f 00 = f°° 
Jo jo 


Thus, with j = 0 and identifying -a 00 /u 0 as angle of attack a, one 
finds from Appendix B 


C La= " 





0 < t 0 



M L* 


cos 


-1 


Mt 0 -1 

■to 


+ — cos” 1 (M 


£%) 


yt 0 2 - ( 1 - Mt 0 ) 



4a 


M+l 


+ 2 ^ r 


(M-l)t 0 






t Q > 


1 

M-l 


Next, with j = 1, and using C^ 1 to designate the pitching moment 
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measured ahout "the leading edge of "the "wing. 




0 ■< t Q 



S{i 




o \ -1 Mt 0 -1 m 

' cos x — - — - 1 - — 


+ — COS -1 (M- P^fco) + 


p 


1+Mt r 


Vto 2 - (l-Mt o y 


_ 1 _ 

6 A 


L MU 


+ 3 t Q - (M- l) t Q 3 ^j- 


— < t Q < — 
Mfl M-l 


^ ” P l 1 " 3 Pa) 


to > 


- M-l 


These expressions agree -with those given hy Miles in reference 2. 

The above results can he used to demonstrate the usefulness of 
equation ( 38 a). Taking j = n = g = 0, Z = 1 in that equation gives 


or, for the present case. 


10 _ -oo „oo 
*00 “ ^00 “ x io 


10 00 00 

f = f - f 
■ L oo x oo 10 


which represents the equality 


V - % + <V 


that is, the lift coefficient for a pitching wing equals the sum of the 
lift and pitching-moment coefficients of a sinking wing (primes indicate 
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the wing is pitching about and moments are measured about the wing lead- 
ing edge). Hence, 



0 < t 0 


< 


1 

Mfl 


2 

M 


1 

K 



COS " 1 


+ M cos-1 ( M - p =t 0 ) 

t 0 e 


+ 


3 -Mt c 


y ± 0 2 - (1 - Mt Q ) : 


1_ 

6a LM+1 


+ 3 t 0 - 3 (M-l)to 2 + 



— < to < -i- 
M+l “ “ M-l 



A further application of equation (38a) provides the pitching- 
moment coefficient for a pitching flat rectangular wing. Thus, with 
I s j s 1, n=g=0, equation (38a) gives 




) 


which becomes 



and so 
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From equation (B2l) in Appendix B it is found that 


,00 

TT 0 ° 

20 

F 20 

^O 

^O 

Uo 

” u 0 


M L3 


( 1 - Mt 0 3 ) 


-^-L- 

12A L 


m(m 2 + 3)t 0 3 


0<t o <1 


= it (if cos-1 !£°± + i M cos -i (M .^ 

M 1« (_ 3 t 0 3 P 


i+Mt 0 + (^+2) t Q g . juf-a. ♦ h. - ‘ 

9 J 24a L I&-1 

(M - 1)3t ° 4 ]} 


Combining, we find 


<V - - - S [ 8 - ® to + M < m2+ 3)to S ]} 0 < to < ^ 

■ - g Jl [ £151°! oos -1 go--. 3 : + 2 M cos -i ( M ) 

M V L 3 to 3 0 

8-Mt 0 -(M 2 +2)t 0 2 / — gl t T „ ^ . n 

— J t 2 _ (1 -Mto) - -f- -Ar + 8t o - 6(M-l)t 0 2 + 

9 J 24A (_ Mfl 

(M - l) 3 t Q 4 1 X — < t Q < — 

= - 2,/i. i_ "l to > -i- 

p\3 4pA J — M-l 


to > -i- 

- M-l 


Another relation among the generalized indicial forces f . can he 

derived by means of the reciprocity relations given in reference 5- The 
details of the derivation are given in Appendix C and there results 


- 1 (S) n 
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Equation ( 39 ) can "be used in two ways; one, as a means for checking 
the internal consistency of a set of calculated generalized indicial 
forces, and the other, as a means for expressing a given force in terms 
of a set of others. 

Consider, as an example of the former use, the case for which 
l = j = 0. Then 


on og 

f = f 
og on 


From equation (37) we can express this relation in terms of the calcu- 

on , 

lated quantities F q ^ thus 



If now n = 1, g = 3 the following relation results 



which provides a useful check on the computed quantities. 


Next let us solve equation (39) tor a given force. Perform the 
sum operation 



on both sides of equation ( 39 ), and reverse the order of summation on 
the left side. There results 
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The inner sum on the left can he evaluated. Thus one has 

P 


- [1 - (l-x)) p = £ (-if (j) (1-xf 


P=° 


-if(?)Z(-« r (0^ r 

r=o 

= £ ( -i)^j(.if( p )( p ) 


Ji=0 

p 


r=o 


Equating coefficients of x, 

P 




ja=r 

and equation (1|0) becomes 


0 r < p 
(-l) P r = p 


& 

Jg 


j=o P=o 


pn 


(41) 


CONCLUDING REMARKS 


A method is presented for evaluating the generalized forces on a 
rectangular wing flying at supersonic speeds and having an aspect ratio 
such that PA > 1. The generalized coordinates used to define the wing’s 
behavior are the amplitudes of downwash distributions expressed in terms 
of polynom ials in x and y, the chordwise and spanwise directions* 
respectively. 

Numerical results are presented in table I for generalized indicial 
forces on a wing having an aspect ratio of k- and flying at a Mach number 
equal to 1.1 and 1.2; the polynomial coverage being 0 < l < 1 and 
0 < n < 5* where w ~ x V 11 * 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field* Calif.* June 30* 1954 
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APPENDIX A 

EXPRESSIONS FOR THE POTENTIAL 


In order to write the expressions for the potential in all regions 
shown in figure 1, it is sufficient to derive in detail only that for 
region V. Having carried out this analysis, one can determine the expres- 
sions for potential in other regions without difficulty. 

Consider, therefore, equation (13) and let a and r apply to region 
V a . First, it is necessary to determine the potentials ¥^ and ¥5 in 
the t,x,| space. From equation (ll), in conjunction with sketch (f ) , 
it is found that 


% = 


. i y Wt2 ' Sl r t ’' /<x " Xl)2+5l£ v u< x i +Mt i’3 r i) 4t i 


X- 


Vt 2 -i x J 


y(t-t 1 ) a -s 1 a -(x-x 1 ) < 


(Al) 


W B = 


n ° t-*/ (x-X-,) 2 * / v 

_ 1 r dx x T v u (x. l + Mt 1 ,y 1 )dt 1 

~ x i/ m y(t-t 1 ) 2 -i 1 2 -(x-x 1 ) j 


*•/ 




\/(t-t 1 ) 2 -! 1 2 -(x-x 1 )‘ 


(A2) 


where 


^1(^1) = ^ ^m ~ J "tm 2 _ ^i 2 ^ 


With the values of ¥ given in equations (Al) and (A2) it is pos- 
sible now to solve equation (6b) for i^, sketch (g) giving the required 
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data in the £,y plane. Thus, if R 2 = (i-t,) 2 - (y-y,)' 


tU>x,y,t) 


at 


Wa 


n J n t TT r>-5+y+t p t 

-*/ **f «nH/ a 

s+y-t s+( y - yi ) y i-(y-n) 



P 2 
-X 


If «*/" 

* i+y-A 2 -^ s+(y-yx) 


-g+yWt 2 -x 2 ,/t 2 ^? 

l-(y-yi) 5+y-t 


i+(y-yi) 0 My+yJ 


a/ t 2 -X £ 


i r ay, r d 5l M + 

l+y-J t 2 -x 2 i+(y-y,) R 


* f 


5-yWi 5 ^ rs^a 


/ 'V O -A 

dl, 

/ \ 


S+Cy+y,) 


W B~ W A 

R 


(A3) 
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Now apply the operation of equation (7) and the potential 9y is 
given hy a 



(A4) 


where Rj 2 = S^ 2 - Cy-y^ 2 and the bars on the integrals signify that 
the finite part of the integral is to he taken in the sense defined 1 in 
reference 10 and that the order of integration c ann ot, in general, he 
reversed. 2 For convenience set 

10 



1 For the subsequent analysis to hold, the definition of the finite 
part given in reference 10 is essential. This definition differs from 
that given hy Hadamard when it applies to multiple integrals. 

2 Since the order of integration plays an important role in the fol- 
lowing development, integration first with respect to x and then with 
respect to y will he denoted J dy / dx f (x,y) while integration first 
with respect to y and then with respect to x will he denoted 
J dx J dy f(x,y). When the notation Jff (x,y)dydx is used, the order of 
integration is immaterial. 
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where In is the nth integral group on the right-hand side of equa- 
tion (Ah ) . 

Consider the first of these integral sets . Using equation (Al) , 
we can write 



In order to simplify this expression, the order of these integrals will 
he rearranged so the integration with respect to can he carried 

out first. The technique of changing the order of repeated integrals 
with strong singularities set forth in reference 10 will he used here. 
Consider the change of order in the x x plane. Pretend for the 

. moment, that the t x integration has 

J i heen carried out. Then the highest 

order singularity (since w u is 
hounded) in the § x , x^^ pla n e has the 
S s ^ order 3/2 which is weak in the sense 

that no residual occurs when the 
/ £/ s y~y/ sequence of integration is reversed. 

Sss/X/ s/ Z /sx/S J The top of sketch (o) shows the area of 

I ■ ’ ► of integration, so immediately 




^ -f dy, / 


x+*/t z -(Y-Y 1 ) i 






M§i 

[S 1 2 -(y-y 1 ) 2 ] 3/2 


t f -t-(x- X t ) 


t-V(x-x ;L ) 2 +i 1 ‘ ! 


w u (x 1 +Mt 1 ,y 1 )dt 1 

' , (t-t 1 ) 2 -(x-x 1 ) 2 -i 1 2 


Sketch (o) 



6M 


MCA TN 3286 


kl 


To change order in the plane, consult the bottom of 

sketch (o). In this case an inherent singularity exists at the con- 
fluence of the singul ar ity lines of the integrand; namely, where 

= y - y x 8X1,1 t x = t - /(x - x x ) 2 + I-l 2 . The change of order can 
therefore not he performed directly, hut account must he taken of the 
existence of a residual term (see ref. 10) . This residual is defined 
as the difference between the two integrals taken in different orders 
over a vanishtingly small region surrounding the inherent singularity 
(the region heavily shaded in bottom of sketch (o) . The residual R^ 
is then. 


Ri = 


lim 

e->0 


/ 


«/(r 0 +e) 2 -(x-x 1 ) ; 


y-y x 


^t-^/(x>j^J 2 ^T 3 2 

/ 'w u (x 1 +Mt 1 ,y 1 )dt 1 

J t-r n -e 


* 1^1 

[lx 2 - (y-y x ) 2 ] 3/2 

t-r Q 


n (t-t-]_ ) 2 - (x-Xi) 2 -!^ 


w u (x 1 +Mt 1 ,y 1 )dt 1 


;-r Q -e 


*/( t-t 1 ) 2 -(x-x 1 ) 2 

/ 


y-yi 


tiUi 

^! 2 -(y-yi) 2 ] 3 / V(t-t 1 ) 2 -(x-x 1 ) 2 - I1 2 




where r 0 2 = (x-x^ 2 + (y-yi) 2 . The second integral vanishes (see 
ref. 10), and, passing to the limit e — >0 in the first integral, there 
results 


jr v u (x 1 +Mt-Mr 0 ,y 1 ) = _ * [w u ] 
2 r Q ~ 2 


where the square brackets again mean that the retarded value is to be 
taken. Thus, the integral I x can be reduced to 




I / W 


xWt 2 -(y-y 1 )^ 


dxn 


y-t 


x-n/tMy-yJ 2 


I>u3 

r 0 


(A6) 
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In the same way, the integral I 2 can he reduced, and 



which is recognized as Kirchhoff’s formula, equation (3).> with an 
acoustic plan form hounded hy the circle 

(x-xj 2 + (y-y x ) 2 = t 2 

The reduction of the integrals I 3 , I 4 , I 5 , and I 6 is quite simi- 
lar, leading to the sum 



(A7) 


Examination of the limits on these integrals shows their total area 
of integration is that shown in sketch (j). But this area corresponds 
exactly to the acoustic plan form S a for a point in region V a l Hence, 

denoting the combination of terms in equation (A7) hy we can write 

simply 



(SaV a 
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b3 


It now remains to calculate the integrals I 7 through I 1Q . Desig- 
nating their total effect on the potential hy one can readily 

show (since no inherent singularities arise in these cases) that 


-y+t x+'A 2 -(y+y 1 ) 2 

^ )= ^ ^ f 


cbc 


r. -1 ” 1 Jbyy x w u (x 1 +Mt 1 ,y 1 )dt ;L 


■f 


[(t-tj.) 2 - r 0 z ]J(t-t x ) 2 - r x 2 


-yW t 2 -x 2 


A/ ^ / 


dx. 


x-^/^-Cy+yi ) 2 


t-r. 


U-J.-L -y+ Vt 2 -x 2 

^ v^yyi w u (x 1 +Mt 1 ,y 1 )dt 1 + ^ f dyi 

o [(t-tj 2 - r^l^t-ti) 2 - r x 2 ^ 0 


t-r 


r’ -1 *Ayy x w u (x 1 +Mt 1 ,y 1 )dt 1 
X 1 (y+y 1 ) - Xl /k [(t-t^-ro^t-ti) 2 -^ 5 


(A9) 


where r^^x-x^ 2 + (y + yi) 2 . Now let 


I 


t-ri 


C(xi,yi) = < 


•J ^yy! w u (x 1 +Mt 1 ,y 1 )dt 1 
~ x iM [(t-t 1 ) 2 -r 0 2 ]V(t-t 1 ) 2 -r 1 2 

^ ri Ayy x w u (x 1 +Mt 1 ,y 1 )dt 1 
° [(t-t 1 ) 2 -r 0 2 ]y(t-t 1 ) 2 -r 1 2 


x, < 0 


X! > 0 


> (A10) 
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In terms of this expression, equation (A9) can he written simply 

4jf = i? f f C(x 1 ,y 1 )dx 1 dy 1 (All) 

( S c ) Va 

where the area (Sc)tt is illustrated in sketch (k) . 

' Si 

In order to give expressions for the potential in every region of 
the wing shown in figure 1, one can show that it is only necessary to 
vary the areas over which the double integration in equations (A8) and 
(All) are carried out. This is evident in connection with the source 
(1) 

portion cp x , for in every case 


cpt 1 ) = - 2. 

2rt 


/ j' 

Sa 


(A12) 


and only the acoustic plan form S a changes with the region. In the 

case of cp( 2 ), the part of the potential due to the existence, of the 
side edge of the wing, equation (All) can be generalized and written 


< T > ^ = -% J'J CU^yJdx-jdy! (A13) 

S c 

where the integrands are defined in every case by equation (AlO) and 
only the "reflected” acoustic plan form S c changes with the region. 

The region Sc is always bounded by portions of the "reflected" circle. 

(x-xj 2 + (y+yj 2 = t 2 


and the "reflected" ellipse 

Xi-xm^ + (y+yj 2 = tm 2 

Figure 2 shows sketches of "both Sc and Sa for all regions in figure 1. 
The absence of a sketch indicates that the corresponding integral does 
not exist for that region. 
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APEEHDIX B 

‘I‘Hh: GENERALIZED INDICIA! FORCES 
The Loading Coefficient 


In order to determine total forces acting on the wing, it is first 
necessary to obtain expressions for the loading coefficient Ap/qQ. 
According to the linear theory 


Ap b 

= UoM St 


(Bl) 


so it is necessary to differentiate each of the expressions for poten- 
tial. As an example, consider, as in Appendix A, just region Va of 
figure 1. The loading coefficient will he divided into two parts 

Ap ( 1 ) /lo and Ap( 2 )/qo correspond to the potentials and 

Thus, using equation (All) 



-y+t x+Jt ^-(y+yj 2 


rt^oM 


/ 


x-Vt^-fy+yJ 2 


§2 dx x 

at 



I 


x “ ''/t^- ( y+y 1 ) ' 


ac 


dx x + 



o 



xjy+yj 


(B2) 


since the derivative passes the x 1 ,y 1 integration without effect. 
Referring to equation (A10) for the function C(x 1 ,y 1 ) we next find its 
derivative -with respect to t. Write t = t-t x ; then for X! < 0 


C(x 1 ,y 1 ) - J' 


VWi w u (x x + Mt - Mr,y 1 ) dT 


(t 2 - r Q 2 ) */t 2 - r x s 
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and 


dc v /, 4yy 1 w^(0,y) 



Notice that if does not depend on (x 2 + Mt x ) the integral term in 

equation (B3) vanishes, -while if it does, then the integrated term is 
zero. Next, for x x > 0, 



In this case, both terms exist unless *w u is not a function of (x 1 + Mt 1 ), 
in -which case the integral vanishes . 

Substitution of equations (B3) and (b 4) into equation (B2) -will no w 
yield an expression for the loading coefficient corresponding to the 
influence of the side edge; 
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^7 


Ap 

9o 


( 2 ) ka 


In 


-y+t 


V fl 


« 2 UoM( 


. Z+n 


x+ 4 t 2 - (y+y ) 2 , 

p dy p Ay y x x^y^dX! + 

o x-«/t 2 -(y+y 1 ) 2 ( t2 “ r o 2 ) A 2 -r x 2 


-y+t x+Vt 2 -(y+y x ) 2 

Ml j dy x j dXj 

0 x-V t 2 -(y+y 1 ) 2 


t-r 


1 V 1 »-yy 1 (x 1 +Mt ;L ) z " ;L y 1 ri dt 1 - y+vx x 


-y+ A 2 -x 2 


^^yy x x^y^dXj, 
. A /t 2 -(y+y 1 ) 2 (t 2 -r 0 2 )A/t 2 -r x 2 


-y+y t 2 -x 2 


- Ml 


dy x 


/ 


t-r 


dx n 


1 a/ 1 tyy 1 (x^Mtj ^) 1 ~' L y 1 n at 1 


-J t 2 -(y+y x ) 2 0 [( t-t x ) 2 -r Q 2 ] a/ ( t-t L ) 2 -r x s 


-y+*/t“-x 


2 „2 


Ml 


/ 


< 3 y 


A 


dX n 


x 1 (y+y 1 ) 


t-r 


P a/ ^yy x (x x +Mt x A 1 y x n dt x 

-(x x /m) E (t - -tL^-ro 2 ]aA^F^? 


(B5) 


The explicit form of given by equation (2) , has “been inserted and 

it is assumed that l > 1. 
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The portion of the loading coefficient corresponding to cpy can he 
found readily and is a 


*MQ 0 c z+n J o 


.y+t (xWt 2 -(y-y 1 ) 2 )V( X -A 2 -(y-y 1 ) 2 ) i 

yi n ■ - 

7t 2 -( y-yj 2 


py+ 1 nx+vf -iy-y x ; 

1 / yi d yi / , 

o x-v / t 2 -(y-y 1 ) 2 


:+Vt 2 -(y-y!) 2 [ Xl +M(t-r 0 ) ] 1 


oy+vt^-x 2 nO 

Ml / y^ay^. / 

0 x-^i^-Cy-y ^ 2 


[ Xl +M(t-r 0 )]' 


dx*. + 


■+ v t 2 -x 2 


p° 

^ / 

x i(y-yi) 


[xj+MCt-ro )] 1 


y+Vt^-x 2 1 

p n (x-yt 2 -(y-y 1 ) ji ) 


'tMy-y^ 


It is clear that, even for small values of the indices l and n, 
the required integrations for the determination of total forces on the 
ving pose formidable problems. There is, however, a property of the 
loading coefficient corresponding to vertical velocity distributions of 
the type chosen here (eq. (2)) that will materially shorten the requisite 
labor. This may be expressed as follows , adopting the convention that 

Ap^ n /qQ corresponds to a downwash distribution proportional to 

(x+Mt^y 11 : 


5 Ap* n 1 Ap Z -^ n 

8x q Q c io 


1 > 0 
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or, 


Ap Zn 



(x!,y,t)dx 


l > 0 


(B8) 


Details of Evaluating -the Generalized Indicial Forces 


In calculating the generalized indicial forces "by means of equa- 
tion ( 36) , it has "been shown that only the value zero need he taken for 
the index l. Thus we must find 


F°“ = - f (x+Mt) J dx f e ay (B 9 ) 

Js w +g+l _ M t i ^ 


The values of the loading coefficient Ap° n /q 0 are found hy differentiat- 
ing the expressions for potential given in the first part of this 
appendix. 

It is convenient, in evaluating equation (B9) , to consider the inte- 
gration with respect to y first. Seating 


L = 




(BIO) 


it is found that L seems to have different representations according 
to the interval in which x lies. These expressions can, however, all 
he expressed hy the same formula. The portions of L corresponding to 

the parts cp^ 1 ' and of the potential are similarly signified, and 

we have 


(1) 




«UoMc 

[n/2] 


n+g 


/ (-if 1 — n l g * — Ko(n+g) + %(n+g) 1 - 

(n+g+l) ! L -* 


EGO 


(s) 


n+g+i -2 p _ 


p=o 


n+g+l-2p L 


Ko(2P-l) + %(2p-l) | 


(2) 


l on 


itUoMc 


J(n^g) 
n+g 2°*® 


^Ko(n+g) + K^j(n+g) 
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■where 


n eos _1 (-x/t) 

Ko(n+g) = t n+g+1 R.P. I sin n+g+1 0d0 


%(n+g) « ^ W 1 * 8 * 1 R.p. J' 


COS _1 (Xjn/tjn) 


sin n+g+1 0d0 


N 11 


. v 2 r 1 dT i r 11 ( T i" T 0 (w,) r-~ w 

J(n,g) = * / 7=% / a - ^ *l x 

J 0 Vl-T| 2 ^ 1-T1. 2 


and [n/2] means the greatest integer contained in n/2. The function 
J(n,g) may he expressed as summations, and it has the property 


J(n,g) = J(g,n) 


(B13) 


The sum formula is, with g + p = n 

[ P/^ 3 / \ r— y . . . . 


i=o 


'p-2i+l 2g+2 \ 


(-D 


g-i 


[p/ 2 ] g-i 

X (Zi)t<-» 3 *( S 1T’*) B ( 

i=o j=o 




p-2 i+2j+3 2g-2j-l 

2 ' 2 


[p/ 2 ] x v i-i 


; I ft) I ■(¥• >0 


p+2g-2i+2j+3 1 


(B14) 
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Values of the function J(g,n) 

\n 

g\ 

0 

1 

2 

3 

4 

5 

0 

m 






D 

1 

^ 3 





2 

5- 8 

— Jt 

4 3 

1 

2 





3 

11 

6 

21 8 

-•7rS+ — 

64 5 

1 

3 

_ 53_ ff . 32 
256 35 




189 w 32 
64 5 

11 

15 

129 128 

256 105 

1 

4 

5329 256 

16384 315 


5 

71 

15 

-J&+& 

256 21 

37 

84 

975 64 

4096 63 

1 

5 

11801 ^ v 512 
" 65536 693 


where 



is the "binomial coefficient 



P* 

(2i)l(p - 2i ) I 


and B(p,q) is the "beta function 


B(p,q) - J (l-x) 9 1 dx 


- 2 f 


It/2 


sin 25 " 1 0 cos 2q_1 0d0 




(B15a) 


r(p) r (q)/r(p + q) 
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The function J(g,n) has heen calculated for g,n taken 0,1,2, 3,4,5. 
Because of the property (B13) , it is only necessary to give a triangu lar 
array, -which appears in the above table. 

Bow consider the functions Kq(v) and K^(v) , defined after equa- 
tion (B12) . It is convenient, for computational purposes, to express 
these in terms of the incomplete* beta functions, defined as 


b 1 _ x 2 (p,q) 



sin 2 ^ 1 9 cos 2q 1 


a-e ) q ' 1 d£ 


0d0 


> (B15b) 


A tabulation of the incomplete beta functions is available in refer- 
ence 11. Note that -when the symbol B is written without a subscript, 
the complete integral is meant, that is, in equation (B15b), x equals 0. 
It is necessary to exercise some care when interpreting Ko(v) and Kj^(v) 
as beta functions because of the upper limit. Thus, since 


, . V+1 ocos" 1 (-x/t) v+ 

Ko(v) = t E.P. / sin , 0d0 


we have the following cases: 

(i) x > t, R.P. cos” 1 = it 


Ko(v) - t W1 B^, 


(ii) 


0 < x < t, R.P. 





Ko(v) = 


2 
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(iii) 


-t < x < 0, R.P. 




lV+ 1 r 


Ko(v) = B 


{ V+2 


2 |_ i-(x/t) 2 ^ 2 * 2 / J 


(iv) -Mt < x < - t; R.P. cos -1 = ° 

Kq(v) = 0 

A similar line taken with %(v) leads to 


(i) x > t, %(v) = 0 

<*, K M (V) = i|tta V+1 


B 


^ L i-(Wtm) 2 


V+2 1 V 

2 * 2) _ 


(iii) - t < x < - %(v) = g f tm +1 


2B 


V+2 
2 * 


£>- 


B 


i-(W'fcm) 2 l 2 ’ 2 


V+2 1 


(iv) -Mt < x < -t, %(v) = M tm v+1 b ( 1 'N 

P V 2 2 J 


The generalized indicial force can now he expressed 


as 


P °n 

Jg 


8a, 


■on /l [ J(g,n) + 2( _ i)n nig! 


jtMU r ,cJ +g+n+1 ^ L 2 g+n 


(n+g+l)! J L 


*Iq (6+n) + 


* Z M (g+n) 


[n/a] 
r~> / 


s g+n+i- 2 |j. r 

1 (*) po (2W - X) + *4 (»-«]} ( B 16 ) 


[1=0 
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where 


_ — XU u p i 

( v ) = J (x+Mt)“ dx t V+1 R.P . J 


cos -1 (-x/t) 


sin v+1 eae (bit) 


*1^ (v) = (x+Mt)^ dx M tm V+:L R.P. J' 


cos 1 (xj a /t m ) 


sin v+1 eae 


(B18) 


It is convenient to express these forces in terms of aimensionless 
quantities. Thus setting 


* _ x + - t 

c' to “ c 


*Iq ( v ) = C J+V+ Y (xq + Mt 0 ) ^ axo t 0 v+1 R.P. 


'cos-^-Xo/to) 


sin v+1 eae = c J+v+2 i 0 j (v) 


*4 (V) = c j+v+2 


i-Mt 0 


f (xo + Mto^dxo mAo^o \ V+1 ep _ 

S V P / 


cos-1 1 

/ *° +Mt ° sm v+1 eae = c 3+v+2 i„J (v) 
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and 
F 


4a,. 


on ‘ Kl on f 1 


J(gjn) 


JtMa 0 l2A L 2g+n 

( - 1)n 2 (^SjT ] [ (e+n) 

[n/ 2 ] /A \ 


g+n- 2 |i 


g+n+l-2}i 


\i=o 


(ii) X 0 = 1 - Mt 0 ; <t Q < -±- 
Mfl ~ ~ M-l 


(iii) X Q = -t 0 ; -L- <t Q 
M-l — 


hMtc 



*7TS2 7T g 

X X 0 - -toCOS# 


( i) Oc 4 <I/(M+I) 


[ I o (2^-D + 4 . } (®2l) 


The integrals Iq(v) and lA(v) 
can "be simplified by reversing the 
order of integration. This can be 
accomplished in a straight-forward 
manner by merely inspecting the region 
of integration in the Xq,6 plane. 
Consider first the integral I^(v) . 
Depending upon the relation between the 
chord length and the time, we see - 
from sketch (p) - that reversing the 
order of integration results in three 
different possibilities for the upper 
limit of the 6 integral. However, if 
we define X Q such that 

(i) X 0 = t 0 ; 0 < t Q < -JL 

~ Lit-1 



/ Mt 0 \ Xo _ _£ C0S Q 

to " • 


-Mto\ 


(it) t/(M-H) <: to <. I/(M-I) 

i X 0 


x a - -A cost? 



hMto -I 




(iii) lAM-l) < t C: 

Sketch (p) 
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then, in 

Iq < V ) = 


every case. 


. V+x r> cos 
/ 

J+l J 


1^ (v) can he written 
~H-*o/*o) 

sin V+1 0d0 


j+V+2 

v o 

J+l 



r=o 



( _ ^oAo) 


sin V+1 0 cos r 0d0 


and, similarly, it can he shown that 


^ (V) 


M 1 
pV+2 j+V+2 



i+MXq / t Q 

M+X0A0 


sin V+1 0d0 + 


Mt 0 J+V+2 

j+V+2 


i 



r=o 



sin 


V+i 


i+i^ M J+i-r 

(B22) 


cos r 0d0 

(B23) 
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APPENDIX C 

DERIVATION OF RECIPROCITY RELATIONS 


According to reference 5, the reciprocity relation for general 
three-dimensional unsteady motion can he ■written 


q 1 (x 1 ,y 1 ,t 1 )W 2 (x 1 ,y 1 ,t 1 ) dx-Ldy^t-L - 
V 

//y> 2 ,y 2 .t 2 )W 1 (x 2 ,y 2 ,t 2 ) dx 2 dy 2 dt 2 ( C1 ) 

V ° 

where the volume of integration V is that swept out in x,y,t space hy 
the wing. The subscript 1 refers to the wing moving in the forward direc- 
tion and subscript 2 refers to the wing moving in the opposite direction 
in the same manner. The coordinate systems are related by 


x 1 = - x e + c - MT 


7 X = ~ y 2 + 2s 


- t 2 + T 


where s,c are wing semispan 
and chord, respectively, and T 
is some fixed value of time. 
These quantities are elucidated 
in sketch (q). 

Now let the wing associ- 
' ated with the subscript 1 have 
the vertical velocity distribu- 
tion 





Sketch (q) 


, j_ \ / X x +Mt i Y/s-yiY 1 

w 1 (x 1 ,y 1 ,t 1 ) =[— c J 


and that associated with the subscript 2 have 


__ ^ N /x 2 +Mt 2 V / s -7z\ s 
w 2 (x 2 ,y 2 ,t 2 ) c J c J 
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Then 


■.WJ • (> - **)'(*?} 

Substitution of these results into equation (Cl) yields 


T c Mti x +Mt 1 \^ o 28 / y x - s \< 

/«*/ *h 1 • -v-\) / ^ C-v-; 

0 -Mt 1 0 


,e A P Zn 

*0 


X C-Mts 7 2S ] 

/»■/ (‘ - a ^) / *. OA) 


Ap^S 


-Mt 2 (C2) 

Equation (C2) can he differentiated with respect to T, yielding 




Ap ln 


c-MT ✓ x 2 + MT\ l P S /y 2 -s V 

-MT o 


y 2 - s Y Ap^s 
<3o 


The binomial expansion is now performed: 
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In equation (C3) the spanwise integration is carried over the -whole wing, 
hut it can easily he reduced to integration over, say, the left panel hy 

use of the factor [1+ (-l)® +n ]/2. Thus, equation (C3) can he written 


(-D 


g 




[l+(-l) g+n ]/2 


c-MT 


sc 


|i=0 


n /x,+MT\H 

I 

-MT 



l 







c-MT 




Ap^® 


By comparison with equations (36) and (37), it is seen that the integral 
terms in the last equation correspond to the generalized indicial forces 

find f , so that the summations can he written 



|i=o 


■I 

1 r=r 



m 


\fheve the quantity (g+n) must "be an even number. 
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TABLE I.- VALUES OF GENERALIZED INDICIAL FORCES, F 

is 

7 n 

The generalized indicial force coefficient F^ is defined hy 

f 

equation (36). It is the response for a mode shape having a unit ampli- 
tude 


a mode = 


x + 


Mty 



and a loading induced hy a unit value of v/u 0 . 


The table gives values of 
traveled) for 



against time (actually chord lengths 


1=0 
3 ~ 0 , 1,2 
n = 0,l,Q,3,k,5 

s = 0,1,2, 3, it, 5 

M = 1.1, 1.2 
A = h 



TABLE I.- VALUES OF GENERALIZED INDICIA!, FORCES, F^ n - Continued 

(a) l = 0; j » 0; M = 1.1 06 



£ 

£ 




TABLE I.- VALDES OP GENERALIZED ODICIAL FORCES, F,_ - Continued 

(l>) l=0; 3 = l; M = l.i jS 



CTN 

U> 


I 
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TABLE I.- VALUES 


OF GENERALIZED INDICIAL FORCES., FrJ 

(c) Z - 0; j - 2; M = 1.1 JS 


- Continued 




£ 




s 


TABLE I.- VALUES OF GENERALIZED IKDICIAL FORCES. F. - Continued 

(d) 1=1; j = 1; M = 1.1 ' 06 



IS 

IB 

1 

£ 

3 

IB 


S 

B 

IB 



4 

5 


IB 

B 

2 

3 

4 

5 

0 

0 

1.212 


1.616 

£.424 

3.879 

6.465 

n 

Iff 

1.616 

£.424 

3.879 

6.465 

n.oe 

2 

IHRH 


3-879 

6.1(63 

n.oe 

19.39 

.055 


1J56 


1.610 

2.429 

3.890 

6.493 

mi 

ippT 

1.61 6 

£.426 

3.883 

6.482 

11.13 


1.6l6 

£.424 

3*801 

6.474 

n.n 

19*4T 

.H 


1UB5 

1.213 

1.622 

2.442 

3.984 

6.575 



1.617 

2.431 


6.532 

11.2 6 


1.616 

2.426 

3-888 

6.530 

n.19 

19-69 

.22 


1.167 

1.219 

1.641 

2.490 

4.044 

6.666 

■ 


1.623 

2.453 

3.970 

6.716 

11.73 


1.603 

a. *i35 

■ESI 

6.604 

EE3 

20.47 

*33 


EE21 

1.232 

1.673 


4-226 

7.290 



1.640 

2.496 

4.080 

6.994 

12.41 


1.639 

2.460 


6-776 

n.94 

21.62 

.ss 


1.169 

1.253 

1.724 

2.675 

4.461 

7.621 


JRv 

1.673 

2.564 


7.338 

13.27 


1.671 

2.510 

4.089 

7*022 

12.54 

£3.07 

.5®4 


1.18Q 

1.209 

1*779 

maul 

4.681 

0.302 



1*713 

2.636 

4.390 

7.700 

14.05 


1*709 

2.560 

4.203 

7*270 

13.10 

£4.39 

*3? 


FW-SU 

HRH 

1.021 


4.861 

8.751 



1.749 

£.701 

4.3E7 

9.013 

14.70 


1.740 

2,616 

V-299 

7.408 

13.61 

£5*6l 

*766 


1.27S 

1.415 

1*997 

3.201 

3.962 

10.24 

* 

5s? 

1.078 


3*017 

9.060 

17-20 


1.069 

2.015 

k.674 

0,276 

15.3s 


l t n 


l t W 

1 

0.1 ' T * 


£ 

11 « 



o.m Nj. 

1 .TOO 

w well 

1ft 05 

IQ ftW 


1 ftftn 

e .mo 

^-063 

rt t n_£ 

17.57 


i.m 


1.330 

1.770 

£•621 


8.2§7 

16.53 


1.691 

2.339 

3^o 

6.861 

13.29 

27*27 


£.306 

3.5J3 

6.021 

11.23 

22.29 

46.59 

2.2 


1.685 

e.000 

3.058 

5*391 

10.49 

21.?7 


1.860 

£.o>6 

4.409 

0.225 

16.63 

35*83 


8.378 

3.9^5 

6-939 

13.39 

£7*71 

60.02 

2*75 


1*796 

2.169 

3.303 

6.m 

12-23 

26.40 


£.017 

2^55 

4.853 

9*e64 

19. e5 

4e.06 


2.773 

4.270 

7-655 

15.03 

31.96 

72. 4l 

3.67 


1.9k 

2.401 

3*052 

7.193 

14.95 

33.77 


2*196 

3.153 

5.492 

10.81 

23.32 

54.14 


3*013 

4.712 

8.599 

17.46 

38-51 

90.95 

5.5 


a.!37 

£.721 

4.519 

0.784 

19.09 

45*27 


a.433 

3.564 

6.395 

13.07 

29.(5 

71*75 


3-376 

5*321 

9*977 

21.00 

S.34 

119*0 

T.333 


e.iU 

2.397 

4.096 

9.738 

21.81 

53.63 


£.561 

3*T9l 

6.90? 

K.4l 

33.39 

04.30 


3*552 

5.657 

10.75 

23*0? 

54.60 

i4o.o 

n.o 


2.31V 

2.996 

5.066 

10.05 

22. £9 

33.9* 


2.650 

3.919 

7.144 

14.ee 

34.05 

85.ll 


3.603 

5*0(7 

n.i£ 

23.84 

56.01 

i4i.7 

0 

3 

2. *24 

3.679 

6.4« 

11.08 

!9*3? 

34.48 

4 

3**T9 

6.465 

11.00 

19.39 

34.48 

62.06 

5 

6.464 

n.oe 

19.39 

34.48 

62.06 

1 12.0 

.05? 


e.4e4 

3^79 

6.W8 

11.10 

19.(4 

34.61 


3*07? 

i 6.465 

11.09 

19.4s 

34.56 

62.29 


6.465 

n.oe 

19.40 

34.5E 

6C.20 

113.2 

.11 


2.426 

3.301 

6.4T9 

n.i4 

19.57 

34.98 


3.081 

6.463 

11.11 

19.4? 

34.79 

62.95 


6.469 

n.09 

19- *3 

34.64 

6C.61 

n4.4 

.22 


sab 

3.996 

6.530 

11.31 

20.07 

36.33 


3*096 

6.493 

U.19 

19*70 

35.63 

65*31 


6.493 

11.13 

19*56 

35*14 

64. n 

1 10.6 

.33 


£.460 

3.936 

6.63a 

11.59 

30.03 

38.31 


3.936 

6.561 

11.36 

20.26 

36.97 

S.70 


6.561 

11. £5 

19*86 

35.97 

66.35 

124.0 

.44 


2*309 

4.015 

6.0OC 

12.00 

21.05 

4o.8l 


4.015 

6.692 

H.65 

£0.96 

38.73 

73.ao 


6.692 

n.47 

20.37 

37.20 

69.55 

132*7 

.jeH 


E.567 

4.109 

6.969 

12.41 

22.01 

43*10 


4.109 

6.849 

11*97 

21.67 

fo.u. 

77.24 


6.049 

11.74 

20. 92 

39.45 

72.53 

140 A 

-?79 


2.615 

4.106 

7.147 

12.78 

23.68 

45*22 


4.105 

6.976 

12. £3 

22.29 

41.93 

00.96 


6.976 

11.96 

21,38 

39*54 

75*23 

146.7 

7w: 


e JVn 


7.7fc» 

Ik. in 

26,69 

52,30 


h.RDl 

T.Wk 


e4,57 

kiMfl 

oq-kfl 


7\«rvk 

12,06 

23,20 


54,56 

i6o*n 

1.0 


3.012 

4.026 

8.399 

15.49 

b - 9 i 

60.03 


4.625 

8.046 

14.36 

06.96 

5e^60 

107^ 


eioU 

13^79 

25^00 

47^76 

94*52 

193*9 

1.371 


3.484 

5.602 

9.?63 

19.02 

38.44 

81.44 


3.359 

9.335 

17*01 

33.04 

67*63 

144.8 


9.S£6 

16.00 

£9.68 

58.42 

120.8 

260.7 

e.e 


3.906 

6.30T 

U.46 

££.62 

(7.61 

105*9 


6.876 

IO.51 

19.54 

39-18 

63*54 

1^7*5 


10.40 

18.01 

34.00 

69.1? 

140.9 

336*8 

2.73 


4.209 

6.ae5 

12.58 

25*34 

54.77 

125*6 


6.771 

IX* 37 

£1.43 

43.83 

95*92 

222.1 


n.32 

19.49 

37*35 

77.34 

1T0.7 


3.667 


4.612 

7.929 

l4.l6 

29*36 

65*76 

177*2 


7.431 

12.54 

24.10 

50.7a 

114.9 

277.1 


1£.43 

21.49 

41.96 

89.34 

£04.2 

405.9 

5-5 


3.130 

8.496 

16. 40 

35.21 

02.24 

206.1 


8.312 

14.15 

27.07 

60.69 

143.3 

362.3 


13*92 

£4.25 

46.47 

106.7 

254.1 

646.7 

7.333 


5*415 

9.034 

!7*66 

36.64 

92*66 

240.2 


6.730 

15.04 

29*99 

66.34 

I6l.fi 

4fil.fi 


14.60 

£5*77 

56,14 

116.9 

£65.5 

750.7 

11.0 


5.626 

9.33? 

18.26 

39-91 

95*12 

243.4 


9.089 

15*54 

31.01 

68.73 

164.1 

427.2 


15*24 

26.63 

53-92 

120.6 

£93-3 

7*1.9 


i i 


I 
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TABLE I.- VALUES OF GENERALIZED JUDICIAL FORCES, F^ n - Continued 

(e) Z = 0; J = 0; M = 1.2 Js 


o\ 

C7\ 


V 

a 

0 

1 

2 

3 

4 

5 


0 

1 

2 

3 

4 

5 


0 

1 

2 

3 

4 

5 

0 

0 

3.333 

3.333 

4.444 

6.667 

10.67 

17.78 

1 

3.333 

4.444 

6.667 

IO.67 

17-78 

30.48 

2 

4.444 

6.667 

IO.67 

17.78 

30.40 

53*33 

.06 


3*293 

3.334 

4.448 

6.679 

10,70 

17. & 


3.333 

V.<*V 

6,671 

10.68 

17.63 

30.60 


4.444 

6.66 7 

10.67 

17-00 

30-55 

53.55 

.1£ 


3.233 

3.33^ 

4.460 

6.713 

10.79 

18,09 


3.330 

4.444 

6.682 

10.73 

17.96 

30.97 


4.444 

6 , 6 ej 

10.69 

17.07 

30.77 

54,15 

.24 


3-J5K 

3.337 

4.499 

6.835 

11.12 

18.91 


3-333 

4.445 

L ».Lm 

6.723 

£ fHh 

10.69 

10.45 

32.2ft 


4.443 

L L Li 

6.667 
<** /'/'-» 

10.74 

-1 A ttA 

18.11 
.a l 

31.55 

56.JV 

:Sa 


3.060 

3*3^ 
3. 3*5 

4.619 

7-wwy 

7.211 

12.14 

£1.50 


3.299 

<i!W6 

£6*7 

no? 

i9.?7 

S^.-M 

36.41 


*.<•37 

6.667 

10.91 

18.87 

>i.Of 

33.98 

39*^9 

63.20 

.545 


3.058 

3.3*7 

4.654 

7.325 

12.45 

£2,51 


3.292 

4.446 

6.805 

11.54 

20.44 

37-6? 


*.*35 

6.667 

10.96 

19.09 

34.72. 

63.3a 

.6 


3.061 

3.369 

4.710 

7.460 

12.77 

23.10 


3.307 

4.474 

6.956 

11.73 

20.93 

36.94 


4.459 

6.T07 

U.06 

1 9-3? 

35.5a 

67. £ 

.8 


3.1V7 

3.52V 

5.013 

8.107 

14.23 

£6.49 


3.*b 

4.672 

7.p6 

12.67 

23.17 

44.39 


<*.6W 

7.002 

11.69 

20, §7 

39.18 

T6.3V 

1.0 


3.266 

3.706 

5.355 

a. 621 

15.83 

30.22 


3.?o3 

V.907 

1*83* 

13.72 

23.63 

50.36 


4.867 

7.352 

12.41 

22.53 

43.20 

66.54 

1-3 


3.53a 

4.153 

6.202 

10.62 

19.94 

6. 08 


3.952 

5.484 

8.993 

16.33 

31-92 

66.04 


5.401 

3.211 

14.19 

26. 6B 

93-V6 

112.7 

2.0 


3.806 

4.551 

6.975 

12.31 

23.9V 

?o.ol 


4.274 

5.997 

10.05 

lfl.78 

J 7 *?! 

61.67 


5.867 

8.97V 

15.81 

30.55 

63-29 

138.7 

2.4 


3.968 

4.833 

7*531 

13.55 

26.92 

57.60 


4,490 

6.360 

10.61 

£0.56 

fc.46 

93.56 


6.192 

9*513 

16.97 

33.37 

70.58 

158.3 

3,0 


4.221 

9.200 

8. £60 

13.19 

30-91 

67.05 


4.738 

6.033 

11.60 

22.92 

48.40 

109.6 


6.612 

10.21 

10.49 

17.00 

00.29 

104.0 

4.0 


<1.533 

5.681 

9.196 

17.28 

39.95 

80.75 


5.17a 

7-V53 

13.08 

25.93 

56.07 

1*9-7 


7.162 

11.13 

20.46 

41.03 

92.59 

210.1 

6.0 


4.6?; 

6.173 

10.06 

18.99 

39-72 

89.60 


5.602 

8.O94 

14.28 

28.46 

61.07 

143.8 


7.771 

12.09 

22.31 

45.08 

102.1 

241.7 

0 

3 

6.667 

10.67 

17.78 

>3.46 

53.33 

94,81 

4 

10,67 

17.78 

30.48 

53.33 

9V.81 

170.7 

5 

17.73 

30.48 

! 53*33 

9V.81 

170. T 

310.3 

.06 


6.667 

10.67 

17.79 

30.31 

53.46 

95.18 


IO.67 

17-78 

30.49 

53. <10 

93.03 

171.3 


17.78 

30.40 

> 53.35 

9V.92 

171.1 

311.V 

.12 


6-OT 

10.67 

17. 01 

30.62 

53.82 

96.22 


IO.67 

17.78 

30.53 

53.58 

93.66 

173.1 


IT .78 

30.46 

53.4a 

95-av 

172.1 

31V.7 

.24 


6.667 

10.67 

17.89 

31.01 

55.13 

99.?6 


10.67 

17.78 

30.66 

5..S3 

97-89 

179*7 


17.78 

30.48 

53.63 

96.35 

176.1 

326.3 

-iC 


6.666 

10* 

IQ rt*l 

31*57 

WT /v\ 

105.4 


in j in 

VT *Tft 

Vi CL 

M 1i 

■k ni i 

■1 Pn 5 


1 T *?.Q 

on kfl 

ni 

<vt rtL 

liH a 

1 

!<5 


6‘.6fo 

10.67 

Jfllii 

32. 2£ 

59*18 

111.7 


10.67 

17*78 

£!cS 

56*23 

104.6 

200.2 


17*78 

Site 

&30 

99 i?9 

i8a!c 

362*8 

.545 


6,664 

10.67 

18. e2 

32*58 

60. 4i 

U5.3 


10,67 

17.78 

131.18 

56.04 

106.9 

206.3 


17.78 

30.40 

9V.50 

100. a 

191.9 

374.0 

.6 


6.703 

10.73 

16.39 

33.06 

61.75 

118.9 


10.73 

!7.09 

! 31.46 

57-66 

109.3 

212.9 


17.68 

30.66 

54.99 

102.3 

196.0 

305*3 

.8 


6.993 

11.20 

19. 4o 

3503 

67.95 

134.7 


11.20 

18.67 

33.18 

61.09 

120.0 

240.5 


10.67 

32.01 

57.97 

109.7 

215.0 

434.7 

1.0 


7*335 

11.76 

e0.58 

38.29 

74.77 

151.9 


11.75 

19.60 

35.16 

66.62 

131.9 

270*8 


19.60 

33.60 

61. Vl 

117*9 

236.0 

488.0 

1*5 


Q.lfe 

13.13 

23.49 

45.19 

92.12 

1 96-9 


13.10 

21,00 

4o.ot 

78.44 

16c. 0 

349.8 


21.6 6 

37.51 

69.91 

130.6 

e09.2 

629.6 

e.o 


8.85© 

14.35 

26.11 

51.61 

100.7 

241. 4 


14.29 

23.91 

; 44.52 

69.43 

190.0 

427.7 


£3.86 

40.90 

77. 6l 

157.9 

34o.o 

768.3 

2.4 


9*399 

15.21 

28.00 

56.28 

121.0 

275.0 


15.12 

£5.34 

47.70 

97.4a 

£12.1 

486.5 


25.26 

43.43 

83. 1J 

171.8 

377*5 

0TB.9 

3.0 


10.05 

16.32 

30.47 

62.43 

137.4 

320.2 


16.19 

£7*20 

j 51.67 

107.? 

24o.3 

563.6 


27-07 

46.61 

90-3V 

190.2 

427.4 

1013. 

4.0 


10.92 

17-79 

33.67 

70.31 

158.1 

377.2 


17.60 

29*63 

57*27 

121.4 

£76.2 

66?. 1 


29.45 

50.79 

99.70 

213.6 

490.7 

U90. 

6.0 


11.8? 

19.32 

36.71 

77.07 

174.3 

417-8 


19.10 

32.18 

62.43 

133.0 

304.3 

736.3 


31.97 

55.15 

100. T 

234.2 

540.5 

1318. 
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TABLE I.- VALLES OP GENERALIZED UTOICIAL FORGES, F^ 11 

(f) 1 - 0; j = 1; M = 1.2 ' Jg 


Continued 


■ 


B 

H 

IB 

3 

4 

5 

IS 

0 

IB 

£ 

3 

4 

3 

m 

0 

1 

e 

3 

4 

5 

0 

♦Go 

•IS 

.24 

:U 

.545 

.6 

.8 

1.0 

1.5 

2.0 
£.4 

?5 

6.0 

0 

Hff 

X.O#H 

1.617 

1*337 

1.312 
1.294 

1.313 
1.403 

1.514 

1.&4 

i.66? 

1.665 

1*659 

1*635 

l*5?6 

1.540 

1*503 

1.489 

1.510 

1.568 

1-756 

1-956 

2.119 

2.3% 

2.697 

3.122 

2.222 
2 . w".-v 

2.E19 

2.209 

2.169 

£.153 

2.124 

2.123 

2.214 

2.357 

2.773 
3.220 
3*566 
4.062 

4.773 
5*510 

3.333 

m 

3*363 

3*389 

3*405 

3*403 

3*434 

i:S? 

5*048 

6.098 

6.925 

8.095 

9*731 

11*24 

5-3J3 

5.343 

5.373 

5.484 

5*645 

5*821 

5*904 

6.026 

6*732 

7.600 

IO.07 

12.68 

14.78 

17*74 

21.00 

25*14 

6.889 

6.919 

9.010 

9.357 

9*831 

10.50 

10.03 

U..19 

13.02 

15.20 

fti.46 

28.26 

52.24 

60.14 

1 

1.667 

1.664 

1.657 

1.627 

1.579 

1.502 

1.469 

1.450 

1.451 
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TABLE I.- VALDES OF GENERALIZED INDICIA! FORCES, - Concluded 

(h) l = 1; j = 1; H = 1.2 06 
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